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Abstract. We reformulate the Casimir force between bodies in non-trivial 
background media. The force may be written in terms of loop variables, the loop 
being a curve around the scattering sites. A natural path ordering of exponentials 
takes place when a particular representation of the scattering centres is given. The 
basic object to be evaluated is a reduced (or abbreviated) classical pseudo-action that 
' can be operator valued, and can be obtained from a classical path integral description. 
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In a recent paper [T], Casimir forces were calculated for a collection of dielectric spheres 



m 
rn 

, in a constant background dielectric media. By evaluating the stress tensor compatible 
^ . with the Lorentz force law [21 [3| , it was possible to calculate the force for the case of two 
O ■ cind three sphere systems as a multiple-scattering expansion. This approach is indeed 
very similar to the scattering methods used recently to evaluate functional determinants 
^ . which give the interaction energies amongst bodies [H El El [TJ IH]. Certain aspects of 
the scattering approach and the use of functional determinants were also present in the 
earlier work [9l [TOl [HI [12]. The ground state energy of the system is projected out from 
an effective action which can be written as a functional integral over the physical fields. 
For an overview of these methods see the recent book [13j. Both approaches share 
in common the formulation of translation matrices that map eigenfunctions between 
scattering centres and the corresponding scattering coefficients for each sphere. The 
approach in [1] differs from the functional determinant approach principally because 
it uses macroscopic quantum electrodynamics [3] to formulate the problem. In this 
framework, quantisation in real media is explicitly performed. Using this method one 
evaluates a force directly from a stress tensor; and it uses the fluctuation-dissipation 
theorem to evaluate the two-point function of the noise currents (which in turn are used 
to evaluate the bilinear correlation functions of the physical electromagnetic fields [5]). 
All of the quantum mechanical content of the problem resides in this two-point current 
correlator giving it as an order h expression. One calculates the physical observable (i.e. 
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the force) which must be free of divergencies from the outset. Whilst the two approaches 
will agree for bodies in vacua to leading order, this will not be the case for bodies in 
real magneto-dielectric media. 

In this paper we study some properties of the force on a general submanifold by 
reformulating the expansion in terms of a set of loop variables. The observation in [1] , 
that the total path length of the particular set of scattering events plays a key role in 
evaluating the frequency integral, leads naturally to the concept of the loop integral of all 
the scattering events. This is because in the last stage of the calculation when evaluating 
the force on a body or at a point, we must set the arguments of the Green's function to 
coincide, thus closing the path to form a loop. It can be thought of heuristically as a 
causally propagating photon scattering off the spheres and finally arriving back where 
it started. In a similar way, this is basically what is done when one attempts to evaluate 
the density of states. Loop variable formulations are important in other areas of physics, 
notably lattice gauge theory where the Wilson loop \1A\ plays a ubiquitous role, and loop 
quantum gravity [I5], where the definition of a loop derivative also appears. 

It is possible to further rewrite this simply in terms of a loop integral of the 
conjugate generator of translations and the potential that captures the system (e.g. step 
functions representing spheres, smooth short distance or singular type potentials). Note 
that this is not the momentum, but rather the generator of translations (with dimension 
inverse length), a purely classical operator resembling a covariant derivative and 
therefore without Planck's constant or a mass scale entering at all. Loop integrals have 
appeared before in [T6l [T7] where Casimir energies on different spaces are approximated 
from a semi-classical description (for a description of the density of states in terms 
of Maslov indices see [iH])- This is, practically speaking, equivalent to the optical 
approximation [T9| [20| |2T] in evaluating the Casimir energy and the inclusion of closed 
orbits as opposed to only the periodic orbits in [161 [17]. These works share in common a 
saddle point approximation of some functional integral. Complimentary ideas are found 
in the worldline formalism [22] applied to scalar fields that enables Casimir energies to 
be calculated numerically. For applications of these methods see [20l EU [221 l23l [211 [25] . 

Looking further ahead, we might expect this potential is related to some background 
3-metric or affine connection that depends on frequency, since the response functions are 
non-local in time but may be local in space. Thus an expression similar to the adiabatic 
invariants and WKB/Sommerfeld quantisation may be anticipated, except that this is 
purely classical and contains the length of the loop in an implicit manner. For a related 
earlier discussion of such integrals that are not periodic and turn stochastic see |26j . 

In principle the idea presented here is an exact method; one uses basically the 
Fradkin representation [271 [211 [29] of the Helmholtz operator for the physical fields 
which can then be turned into a particle type path integral. To then evaluate this 
requires either very simple media or approximate techniques. They are expected to 
coincide with those found in [161 [13 [IS] when they are approximated, for example by 
straight line rays, and one could imagine that they provide a useful 'discrete' way of 
realising them. 
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The outline of this paper is as follows. In Section |2] we give a description of the 
field equations we are considering and their stress tensors. In Section [3] we show how 
the A^-sphere Casimir force can be reformulated as a loop integral, together with some 
details of the classical generator of translations and the use of representations. We 
also discuss the overlap with a path integral representation. In Section H] we give two 
examples; electromagnetic fields in a magneto-dielectric background in three dimensions 
and a scalar field in two dimensions. Finally in Section |5] we conclude. 



2. Some generalities 
2.1. Equations of motion 

We consider some general background media described by a set of potentials, through 
which a particular spin state propagates. The three most usual are spin 0, 1/2, 1 (scalar, 
fermion, vector), with equations of motion given by 

{V' + e)<p = (1) 

iYid, + h)i; = ji,'/'\ (2) 

V A /i-^ ■ V A E - (a;Vc^)e • E = iwjSi^ (3) 

V A E + iuB = 0. (4) 

To evaluate the force on some object (sub-manifold O) we need to construct a stress 
tensor out of the physical fields and integrate this over the the volume of the submanifold 

F'{0) = [ <rx{\/jT^'). (5) 
Jo 

If the object is compact then we can calculate a force as opposed to a force density when 
one or more dimensions of the subspace are non-compact. For the three most relevant 
cases, the stress tensors one can consider are 

t;-"'- = 5,09,0, (6) 



^sp^or ^ ^^^^ 



It ^ 

2 

7;r-=E.E,+B.B.-i^«(E.EH-B.B). (8) 

At this point the classical field equations need to be solved subject to the usual continuity 
and boundary conditions. The bilinear physical fields in the above can be written 
in terms of classical Green's functions together with the appropriate noise currents 
that persist in the background. Using the fluctuation-dissipation theorem, the vacuum 
expectation values of the two-point current correlation functions can be evaluated thus 
giving an explicit expression for the stress tensor. 
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3. The loop integral variables 

3.1. A loop integral for the inter- sphere force 

We start with the expression given in [Ij for the Casimir force on a sphere due to 
multiple-scattering for the collection of A^-spheres (all are dielectric). See [1] for a 
related discussion in terms of energy functionals. The physical dispersion force carriers 
are the electric and magnetic fields. For all the diagrams of order N in the number of 
interactions, the A^-body force was found to be (where the frequency is rotated to the 
imaginary axis Q = iu) 



Briefly the terms occurring in Equation (Q are : the difference vectors r[i,j] connecting 
each of the sphere centres; the Mie scattering coefficients of each dielectric sphere 
a^{QR[i]); the translation matrices A^'^ {r[i, j]) that map TE and TM modes from one 
sphere origin to another; the S0{3) vectors |1) that truncate the matrix representations 
to some finite order in the angular momentum L; the derivative operator V[.'j^^*j^j that 
acts on the last translation matrix of the order N polynomial; and the function W 
is a measure due to the two separate contributions in the stress tensor, arising from 
the individual field components EmE„ + BmB„ and the part where traces have been 
performed 5mn(|Ep + |Bp). Of course, this could have equally well been done for spin-0 
or spin-1/2 massless fields that are then the effective force carriers. Obviously the total 
force requires adding up all the lower order disconnected diagrams as well. 

It can be seen from Equation ([9]) that the central object of any Casimir interaction 
calculation is the evaluation of the quantity Z[a, A, W]. In fact as will be shown below, 
this is no more than a specific representation of an object that shares some structure 
with semi-classical physics. It is clearly some type of loop integral, as we start and finish 
on the same point of the particular sphere where the force is being evaluated. We can 
write this in a shorthand notation by the introduction of a path ordered exponential 
(see |14j for details of how this is used in gauge theory) and replacing the truncation 
vectors |1) with a trace over S0{3) indices (up to some integer Lmax), together with 
retaining only the polynomials in the translation matrices that form the closed loop 
index wise 




J2 A''\v[l, t])a'{nR[{\)A^'^ir[t,j])a^{nR[j]) ■ ■ ■ VrKi]A"^'^(r[m, 1])] 



i,j,---,m=2 



ikR[l])-j{R[l])-n{R[l])Wm[l])\l) 




/ N N 



) 



-1 loop{A) 



Z[a, A,W] = Tr Vexpij^Yl a\^R[^)A^''{r[3, i]) 



(10) 
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This means that the first contribution retained is the second order contribution in 
different scattering coefficients (the famihar two body forces). The necessity of path 
ordering in an expression hke this is due to the fact that in general we are dealing with 
non-commuting matrix representations. Thus to form a continuous path we must specify 
the order of the translation matrices so that they are sewn together in the correct way. 
Indeed, Equation f lTOj) is really defined by the power series expansion, as the exponent 
is not meaningful for a set of translation matrices of a different dimensions (for example 
a dipole quadrupole interaction). The suggestive feature of this form however, is that 
we encounter a closed loop of finite distance translations in an exponential similar to 
what is encountered in semi-classical physics. 



3.2. The translation operator 



We now try to develop the loop integral formulation and see the connection with path 
integral methods. Define the translation operator 

ry 

dz'DJz) 



U{y,x;D) := exp 



11 



where Di is the classical generator of translations in a given background. The translation 
matrices that map the vector wave-functions between different scattering centres [30] 
are then just representations of this operator in the 5*0(3) TE/TM basis. For example, 
the A matrices used in [H |30] are just the constant finite displacement versions of this 
operator (the background is constant, the spheres are viewed as perturbing potentials, 
and Di = Vj is the standard derivative operator) 



AilZlik\y-x\) 



l^,{L2,m2,hl^{k\x\)\U 


■U{y,x;D)- 


u 


Ll,ml,hl,{k\x\)) 


[J^,\{L2,m2,hl,{k\^\)\U 




u 


Ll,ml,ht,{k\m'V^' 



.(12) 



The expression we have encountered for the force in Equation ffTOj) can be partially 
recast using Equation (fT2l) . By using the orthogonality of the vector eigenf unctions, we 
see that all we are doing is furnishing a representation of the translation operator by 
the insertion of complete sets of states when we have perturbing scatterers present. 

More generally, if we were now to change from the step function potentials relevant 
to the hard sphere scattering (together with the application of continuity equations 
across the surfaces) to a general 'softer' set of potentials and consider Z at a particular 
point (or on a submanifold) , then we would not be dealing with straight line rays and 
reflection coefficients, but rather curved rays in some effective potential. This is a further 
prompting to write the fundamental objects in terms of loop variables. For a spin-s field 
we write the corresponding loop integral as 



Z^'^[n,x] = (F^lexp 
= Tr ( exp 



dq'D^{q,n) 
-Sr{Cx, fi) 



{I,F,x,s} 



(13) 
(14) 
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where x is the initial and final point of the loop, and |I) and |F) are the initial and 
final states (which may differ for example by the boundary conditions imposed on the 
eigenf unctions used to represent in and out modes). The generator Di of translations 
and the loop C^^ implicitly depend on the background potentials. In the next subsection 
we shall show how this loop integral arises in a path integral description. The canonical 
generator of translations Vj (i.e. the gradient operator for smooth functions) is promoted 
to a constrained covariant derivative operator in general 



a dimensionless object, since it is defined using the generator of translations rather 
than the canonical momentum. It is in general operator valued, being realised on 
eigenf unctions of Helmholtz type operators. Obviously it is a functional of both the 
background effective potential and the loop. This should be contrasted with the semi- 
classical evaluation used in [161 [H] where a saddle point approximation is used to 
evaluate the density of states. It is also different from the optical approximation |19j . 
where straight line elements dominate the path integral (giving a stationary phase) of 
a classical particle action with the mass scale set equal to unity. Since the loop integral 
developed above uses the fluctuation dissipation theorem to evaluate the quantum part, 
the physical flelds 'momentum' density will diminish with subsequent windings of the 
loop, thus deflning part of the perturbation theory. The tangent ray in contrast does 
not have this dynamical information encoded. Rather than considering periodic orbits 
in the usual particle phase space, one would expect the conflguration space of windings 
is the appropriate one for dissipative backgrounds. See also [26] for earlier work on 
non-periodic orbits in which trajectories become stochastic. 

The last part to fully specify the Z function is in flnding a set of loops in the 
background that describes how the force carriers propagate. This, together with the 
corresponding representations of the translation operator (if it necessary to cut the loop 
integral up) are the remaining quantities required in any evaluation. The assumption 
we now make in order to calculate the loop equation is that the background effective 
potential can be turned into an affine connection with which to deflne parallel transport. 
Thus the loop is given by the parallel transport of the loops tangent vector based on a 
connection derived from the background 



The significance of this statement is that it is only when the background potential 
possesses some type of critical point, singularity or sharp boundary that the loop integral 
can be non trivial. Equation f|T6|) needs to be solved subject to the boundary conditions 
q^{initial) = q^{final) and q'iinitial) = q^final). It thus tells us the shape of the loop 
rather than its length, since we are not initially assuming a metric connection. In this 
manner once a connection has been derived, one is able to work out the space of loops 




d 

dr 



4 + ^d'q' 



0. 



(16) 
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(the set of closed paths labelled by r G [0, 1]) 

= {q%T)\ if + r;.,g^g'= = 0, g(0) = g(l) = x, g(0) = g(l)}, (17) 

which in turn can be used to then evaluate the loop integral. Note that this prescription 
is different from choosing the initial and final momenta to coincide due to the dissipation 
in the background. Here instead we make only a statement about the shape of the loop 
rather than the dynamics. 

The total loop integral must take into account all possible interactions which is 
reformulated into all possible periodic windings of the loops around the critical points 
(or non trivial features) of the background. One might think that there are a large 
space of such loops with equivalence under smooth deformations (i.e. not crossing 
critical points boundaries, etc.). The use of a connection to give a geodesic equation has 
the effect of collapsing this space down to just the one element for a given frequency. 
The final expression we find for the total Z function is 



Z^'\n,x] = ^(F^lexp - I dq'Di{q,n) 

n L ^ Cx 

= 5^Tr(exp [-SR{C,,T,n) 



\K) 



. - - {IJ,FJ} 

*~'X 

and the sum over loops can be expanded (for example in terms of winding numbers). 
This expression captures the essential elements of the classical scattering of the physical 
fields, together with Equation where the derivative acting on it removes the non 
interacting pieces. 

3. 3. The relation to the Green function 

The loop integral we have derived may be thought simply to be the Green's function 
evaluated at coincident points. This is not quite the case as can be seen by considering 
a path integral representation (see [22]) specifically for the electromagnetic fields. We 
consider now a more general background media given by the dielectric permittivity 
eij{x, fl) and magnetic permeability fiij{x, fl) tensors, i.e. they have no spatial dispersion 
and are just local functions of position. Rewriting the field equation for the electric field 
as 

(^^e-i ■ V A /i-i ■ V A +1-^ E(x, n) = (^^^ e-' ■ jj^^ (x, ^]), (19) 

we see the natural appearance of a dimensionless differential operator on the left. 
This operator has the same eigenfunctions and eigenvalues as the standard Helmholtz 
operator for the electric field. It can be used in the construction of propagators, 
whilst a modified evaluation of the noise currents two point function can be found 
to evaluate the electric field two point function. The corresponding two point function 
can be represented as (promoting the partial derivatives to covariant derivatives with 
connections based on the permittivity and permeability) 
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A{z, y) = {z\ [-cVl^^V, A A \y) 

POO 

Jo 



(20) 



By introducing an integration measure over paths connecting the two point this can be 
written as (where for simphcity we take the scalar version of this Helmholtz operator 
containing the two potentials) 



dr / [dq{t)]^/d^^i^)^J^^)(^^exp 



.(21) 



Considering now the loop integral version z = y for the closed orbits it is obviously 
necessary to compactify the real line r G 5*^ to the reduce the integral to a sum over 
windings, together with a specification of the integration measure. This is similar to 
the path integrals used in [TB] and |22] where a fictitious time is introduced as well as 
mass scale set equal to unity and an appearance of Planck's constant in the particle 
action. As pointed out in [19j, it is misleading to call this a semi-classical evaluation 
because of the absence of the dimensionful Planck's constant. Concurrent with this is 
the absence of a mass or length scale with which to define a dimensionless action. The 
choice we have made in the above can be seen to be dimensionally correct and relevant 
at the local level. It useful to compare this expression with the first order form of the 
quantum path integral for a relativistic particle [29] in a background metric:- 

A{x,x) := {x\ _2„2 1^) 



dr 



<?[0] =<?[!] = 



[dq] \J (\.ei{g) / [dp/h][de 



exp (^i/h j dt[qp + e{p^ + w?(?)]^ . 



(22) 



Performing the integral over momentum, we get back to the configuration space (a 
proper treatment also requires gauge fixing the world line metric e = 1 and the 
introduction of ghosts). It has the feature of inverting the background potentials used to 
define the Hamiltonian (i.e. g^'^p^p^ + m?c^ becomes g^vq^q^ + /rn^c^). If we perform 
first the integral over the world line metric, the mass shell constraint is implemented 
(see [35] for a version implementing Dirichlet boundary conditions in field theory) 



A(x, x) 



dr 



q[0]=g[l]: 



[dq\\/ det{g) j [(ip//i]exp (^i/h j dtqp^ 5{p 



2 + mV)(23) 



(i.e the vanishing of the Hamiltonian on physical states) leaving the form of the reduced 
action operator encountered previously in Equation (fT8|) . Integrating over both the 
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momenta and the world line metric would give us the nonlinear form only depending on 
the position 

A{x,x) = / dr (J) [dq]^/det{g) exp I —irric/h / dt\/ QijCf-q^ j .(24) 

JQ Jq[Q\=q[l\=x \ Jo J 

The significance here is that the exponential contains both the path length of the loop 
and the perturbing potentials via the metric. Expanding the metric about the given 
background case, the perturbing potentials act as the scattering elements. The final 
schematic form then we take for the -E-loop integral is 



Z[VL,x\ = (F^I / dr i [dq]yjdet{g) i [dt)][de] 

Jo Jq[0\=q[l\=x J 

■ exp ^ dt[qb + e{g'W,D, + 1)]^ (25) 

where Equation 0181) can then be considered a shorthand notation for this. The classical 
equations of motion then lead to closed geodesies defining the loops. In addition the 
integration over the frequency can provide the sum over all paths (i.e. a sum over 
all loop sizes where there is one geodesic for one particular frequency), whilst the tau 
integration once compactified gives the windings of the loops. For the scalar version 
given by Equation ( |2T|) this is 

Z[n,x\ = ^(F^l / dr i [rfg]Vdet[fiVc2(e-/i)] / [dV][de] 

n=l -^0 Jq[0]=q[l]=x J 

■ exp (^-n ^ dt[q ■ V + e(V, ■ + 1)]^ (26) 
dVLZ[VL,x\^^Z[VL,x\. (27) 

locyps 

4. Some Examples 

4.1. Electromagnetic fields in media 

In the case of electromagnetic fields propagating in a general background media, the 
loop integral is given by 



and 



Z^[n,x] = (F^lPexp 



(f dq'bj 



%)■ (28) 



To recover the previous results for the A^-sphere scattering in Section [3], one can pass to 
the nonlinear form similar to Equation f l2^ and perform a perturbative analysis. The 
first step would be to expand the material properties around a background value to 
linear order in the fluctuations 

e(g,fi) = (eB(g,^]))+5e(g,fi), (29) 

/i(g,fi) = (/iB(g,fi)) + 5/i(g,fi). (30) 
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The fluctuations then describe the presence of each of the spheres, and the nonhnear form 
can be expanded (the square root of the corresponding hne element). Following this it 
is necessary to introduce complete sets of vector eigenfunctions (spherical harmonics \L) 
and Bessel functions \zi) that solve the Helmholtz equation for the constant background) 
centred on each of the spheres (coordinate system Sj) 

n = {\TEi) = L|L, ZL, Si), \TMi) = V A L|L, zl, S^), 

|L0,) = V|L,Zi,S,)}. (31) 

Both the E and B are expanded in the complete set of states H with no longitudinal 
component to satisfy the divergence free equations of motion. The loop integral can be 
cut up into essentially all the path ways connecting the sphere centres. By inserting 
the complete sets of states we extract both the translation matrices already given in 
Equation ([T2|) and also the scattering coefficients (where k = ^/eBflB^/c) given by 

al = -^ {TMMTM,), (32) 



B2 



t.n2 r 

/5I = ^ JJTM.,MTM,). (33) 

The canonical momentum density (that is consistent with the Lorentz force law) is purely 
in the longitudinal direction V, which consists of the scattered modes with respect to 
sphere-i 

P = E"^ AB, 

= -E"^ A V A E, 

= ^V(Et.E). (34) 

Since it is only the derivative operator entering here giving the direction of the 
momentum, this appears at least to be consistent with the simple derivative operator 
evaluated on the loop. The choice of physical momentum is a subtle point [3ll |32] (see 
also [331 El]) which we do not address in detail here. 

4-2. A scalar field in a potential 

As a toy model consider a pseudo-Hydrogen atom in two dimensions, with a scalar field 
as the force mediator. That is we have two polarisable particles, one of which is big 
(the nucleus) compared to the other (the electron) that are chargeless and interacting 
through the dispersion force carrier (the scalar field). It should be stressed that whilst 
this is not a physical example, it can give some insight into handling singular potentials. 
The .Z-loop integral can be written as 



Z°[n,x] = ^(F^.|Pexp - / dqW 



, l^"^- ^^^^ 
The system has a U{1) symmetric potential /c^(r, Vt) with a length scale R given by 

fc^(r,(]) = ^ = ^m (36) 



Casimir forces from a loop integral formulation 



11 



and is compatible with some two dimensional deformed metric 

ds^ = V{r){dr^ + r^de^). (37) 

The geodesic equations can be simplified using the constants of motion. Due to the 
U{1) symmetry and the invariance of the action in the label r we have 

= 0, (38) 

dr 

-^{V[r^ + r^e^]) = 0. (39) 
dr 

These can be rearranged such that they take the form one usually encounters in classical 
orbits, 



d^u 



EdV 

2 AT- («' 



where u = 1/r and E is the integration constant associated with the constant energy 
orbits. Given the long range potential in Equation f p6|) . it follows the loops are ellipses 
where the eccentricity S is determined by the frequency Q, the length scale R and the 
integration constant E (which can be chosen to be unity). The final sum then needs to 
be over all ellipses with the eccentricity being summed/integrated over, together with 
the winding number n for how many times we wind around the ellipse 



n=l S 

exp[—l{x, £)] 



^ - exp[-l{x,£)y 



(41) 



Here / is the length of the elliptical orbit that depends on the eccentricity and the 
point X where the loop integral is being evaluated at (which determines the integration 
constants in the solution of Equation fHOj) ). 



5. Conclusions 



The principal result of this paper is that we can represent Casimir interactions (which 
are given by bilinear correlation functions of the physical fields) in terms of an equivalent 
particle type path integral with an underlying loop structure . In passing from a field 
description to that of a particle description, we also find that the material properties 
(for example the permittivity and permeability) in which the fields propagate are 
turned into an appropriate geometric quantity (a connection or metric) in which the 
particle like quantity moves. By introducing the use of a connection to define parallel 
transport, it is possible to calculate the space of loops, which can be found from a simple 
geodesic equation. One can pass between different representations of the path integral 
by integrating over the worldline metrics and/or the conjugate translation operator 
(which is locally defined). From here we see how the loops come out from the second 
order form, how the nonlinear form can be used to evaluate expressions perturbatively. 
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and how the first order Hamiltonian form has the translation operator contained within 
it. The final integral over frequency produces a sum over all the loops and in this way 
one finds a sum over all paths. It is expected to be principally useful when applied 
to Casimir physics, for those situations where the materials between bodies have real 
length scales associated with them; at certain frequencies one can imagine being able to 
resolve these length scales and the orbit sees an effective background curvature in which 
it is moving. To go beyond a straight line approximation in this situation is necessary 
and we expect the loop integrals presented here to be applicable. 

It would be interesting in the future to investigate further the role of loop variables 
as the real dynamical variables from a loop quantum electrodynamics perspective [36] . If 
we write the dynamical variables explicitly in terms of the loops, we may well expect all 
the physical observables to be finite and not require the type of regularisation usually 
associated with quantum field theories. The inter body force is one such observable, 
but possibly one could apply it to constructions involving friction [37] where one would 
have to use a theory of dynamic loops. In addition it would be interesting to see how 
the functional methods developed in [3] could be related to the loop integral picture 
presented here. This will be left for future work. 
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